Electron-drift driven ion-acoustic mode in a dusty plasma with collisional effects 
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Instabilities of ion-acoustic waves in a dusty plasma with electron-drift, collisional, and dust 
charge fluctuations effects, have been investigated. The regimes are clearly marked out where the 
theory is applicable. The critical electron-drift velocity required to drive the instability is predicted. 
It is also shown that electron thermal conductivity and charged grains concentration enhance the 
growth of the ion-acoustic mode whereas ion-viscosity, ion-thermal conductivity, and dust charge 
fluctuations have a stabilizing effect. 



I. INTRODUCTION 

Dusty plasmas have acquired considerable importance^ 
because of their applications to astrophysics and plan- 
etary physicsH, problems of plasma processing) 3 ], and 
physics of strongly coupled systemsu. Experimentally, 
one of the important areas of investigation is the study 
of low frequency (ui < fli ) fluctuations in a plasma with 
dust grains. In a recent experiment by Barkan et al.B, it 
is shown that the phase velocity of ion-acoustic wave in- 
creases with the density of negatively charged dust grains 
and thus ion Landau damping hecomes less severe in pa 
plasma with charged dust grainsu or with negative ionsQ. 

The ion-acoustic wave is a typical compressional mode 
in which ions provide the inertia and electrons bring in 
the pressure effects through their shielding cloud. In a 
magnetized plasma, it propagates with acoustic phase 
speed along the magnetic field lines. Earlier, the ki- 
netic theory of the current-dri'umi ion-acoustic instabil- 
ity, in a fully ionized, collisionallffl and collisionlessa two- 
component plasma, has been investigated by various au- 
thors. In the limit of strong collisions (i.e. mean free 
path smaller than the wavelength), Coppi et al.E3 and 
Rognlien et al.Ea have used the two-fluid equations for 
an investigation of the cupent-driven ion-acoustic insta- 
bility. Kulsrud and ShenE3 have investigated the effect 
of weak collisions on ion-acoustic wave using a Fokker- 
Plank model, while StefantEJ used kinetic theory and a 
collision integral of Bhatnagar-Gross-Krook type to show 
that ion-viscosity leads to a damping of the ion-acoustic 
instability, whereas, the electron- ion collisions tend to in- 
crease the growth rate. Recently, RosenbergQ has carried 
out the investigation of current-driven ion wave instabil- 
ity in a fully ionized collisionless dusty plasma using ki- 
netic theory to predict the critical electron-drift velocity 
relative to the ions, required to drive the instability. 

Interpretations of ion-acoustic wave excitation in a 
dusty plasma experiment have typically relied on the 
theories .based on collisionless inverse electron Landau 
dampingt2rliil. A realistic analysis of the experimental 



situation shows that in a typical experimental setupBO, 
the plasma density may be quite high ~ I0 9 -10 10 cm -3 
when the plasma temperature is quite low. As an ex- 
ample, when plasma density ~ 10 10 cm -3 , with elec- 
trons and ions having approximately equal temperatures 
T e ~ Ti ~ 0.1 eV, the electron-ion collision mean free 
path (A e ) may be comparable or even smaller than the 
plasma scale length (L) along the magnetic field, i.e. 
A e < L. Under these conditions, the electrons no longer 
follow the Boltzmann relation and a collisionless theory 
of the electron-drift driven ion-acoustic instability is not 
justified. It is therefore important to re-examine this ion 
wave instability in a collisional dusty plasma. In these 
experiments cited above, however, ion-acoustic waves are 
launched by applying a sinusoidal voltage to a grid im- 
mersed in plasma rather than by drawing an electric 
current. Nevertheless, the effect of Coulomb collisions 
in presence of current driven ion-acoustic waves in such 
a plasma, as is considered in this work, remains to be 
seen in a laboratory experiment. In another experiment 
by Barkan et al.Ej, current driven ion-acoustic waves in 
a dusty plasma are studied. But, this experiment was 
conducted at a relatively low plasma densities ~ 10 6 - 
10 7 cm~ 3 . 

In this paper, we investigate an electron-drift (or elec- 
tron current) driven ion-acoustic mode in the presence of 
negatively charged dust grains. The dust charged grains 
are considered to be massive particles in a multicompo- 
nent plasma. This is valid when a <C d <C Xd, where 
a is the average dust radius, d is the average distance 
between the dust particles, and Xd is the plasma Debye 
length. 

The paper is organized as follows. In Sec. II, the 
regimes of interest, where the theory is applicable are 
discussed. Sec. Ill deals with the linear instability theory. 
Sec. IV contains the conclusions. 
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II. REGIMES OF INTEREST 

At first, we refer to regimes of interest where the 
present calculations are relevant. We consider a typi- 
cal experimental situation-such as considered in exper- 
iments by Barkan et al.B-EJ and Song et al.ti In such 
conditions, the typical parameters in these experiments 
are : potassium plasma (K + ) with mi/m p ~ 40, mj and 
m p being the ion and proton masses respectively, dust 
mass rrid ~ 10~ 12 gm, gas pressure < 10~ 5 Torr, the 
corresponding neutral (n) density n n ~ 2 x 10 cm , 
equal electron (e) and ion (i) temperatures T e ~ Tj ~ 
0.1 eV, and dust grain temperature Td ~ 0.03 eV (as- 
sumed to be at room temperature). The plasma is 
confined radially by a magnetic field B ~ 0.4 T and 
the plasma length along the magnetic field L ~ 40 cm 
with a radius r p ~ 2 cm for the plasma column. The 
characteristic frequency of the waves / ~ 20-80 kHz 
which corresponds to w = 2irf ~ l-5xl0 5 rad/s. The 
other plasma parameters are : electron-thermal veloc- 
ity c e = (Te/me) 1 / 2 ~ 1.3 x 10 7 cm/s, ion-thermal ve- 
locity a = {Ti/nn) 1 / 2 ~5x 10 4 cm/s, the ion-acoustic 
speed c s = (Te/mj) 1 / 2 ~5x 10 4 cm/s, the dust-thermal 
velocity Cd = (Td/md) 1 ^ 2 ~ 2cm/s, the ion gyrofre- 
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the ion Larmor radius 



quency f2j = eB/rriiC 
di = Ci/fti ~ 7 x 10 -2 cm, the collision frequency of 
electrons and ions with the neutrals v en = n n a en c e ~ 
1.3 x 10 4 s" 1 , = n n cr 4 „Ci ~ 50 s^ 1 (where cr in ~ a en ~ 
5 x 10 _15 cm 2 ), with the corresponding mean free paths 



A 



A,- 



en r ~ 'Hn r ~ 10 3 cm. The e-i collision frequency {v e i) 
and i-i collision frequency (vu) in the form presented 
by BraginskiiH are calculated for the present parameters 
and are found to be as v e i = l/r e ~ 8.4 x 10 6 s _1 and 
va — 1/rj ~ 2.8 x 10 4 s~ 1 and the mean free paths are 
A e = c e /v e i ~ 1cm and A; = Ci/vn ~ 1cm, where r e = 

3.5 x 10 4 T e 3 /2 /nA, n = 3.0 x 10 6 (rn,i/2m p )^ 2 jf /2 /nX, 
and A = 23.4 - 1.15 log ra + 3.45 logT - 9.5 (for T < 
50 eV) is the Coulomb logarithm. 

From these calculations the following assumptions can 
now be made 

(1) the collision mean free paths of electrons and ions 
are comparable to or even smaller than the plasma 
scale length along the magnetic field i.e. A e ,i < L. 
so that we can make use of Braginskii's fluid equa- 
tions, 

since u < f e ,-, we can neglect the electron inertia, 
since <Zj < r p and assuming a homogeneous plasma, 
the ion motion transverse to the magnetic field can 
be neglected (the polarization drift), 

(4) since v el > v en ,v ed and i/ u > v in ,v id , so the mo- 

mentum and energy losses of electrons and ions to 
the neutrals and dust grains can also be neglected 
in the fluid equations, 

(5) since the phase velocity of the observed ion-acoustic 

wave is much larger than the dust-thermal speed 
i.e. v p 3> Cd, we can neglect the dust dynamics 
in the analysis. However, we shall incorporate the 
self-consistent fluctuation of charge on dust grains, 



which is shown to .cause a damping effect on the 
ion-acoustic waveHLj. The charge fluctuation on 
the surface of the dust grains, which depends cru- 
cially on the dust size and the plasma density, is 
important when charging rate is comparable to the 
growth of the mode. 
(6) we consider negatively charged grains as an addi- 
tional charged plasma species of uniform massive 
particles similar to a multicomponent plasma. This 
remains valid as long as a <C d <C Xd- 



III. INSTABILITY ANALYSIS 



We write down the basic equations in order to de- 
rive the linear dispersion relation of electron-drift driven 
ion-acoustic wave instability including the effects of 
dust charge fluctuations, electron and ion temperature 
perturbations, ion-viscosity etc. We consider a one- 
dimensional, plane, homogeneous plasma with low j3 
(j3 = 4irnT e /B 2 1 is the ratio of kinetic energy to mag- 
netic energy) . The resultant governing equations for elec- 
trons and ions are the Braginskii's fluid equations viz., 
equations of continuity, momentum, and energy. 

We write down the equilibrium equations that define 
this specific regime of interest and confine ourselves to 
the parallel motions only. They are the ion and electron 
continuity equations 



dt 
dn e 



+ V|| ■ (n<«t||) = 0, 
+ Vii • (n e « e |i) = 0, 



parallel electron and ion momentum equations 

dUi\\ 2 

■ -V||Pi + emE + /i||V||tii||, 

"Viipe — en e E + R, 



dt 



du. 



ell 



dt 



(1) 
(2) 

(3) 
(4) 



where R represents the momentum rgained by the elec- 
trons through collision with the ionstl, 



R — Ru + Rt , 
R u ~ —m e n e v e i(0.51u\\), 
R T ~ -0.71n e VyT e , 

till = U„\\ — Will. 



(5) 
(6) 
(7) 
(8) 



In the above equations, E is the electric field driving 
the instability and flu = 0.96?i;o2^o/j<m is the parallel 
ion- viscosity coefficients. In the ion momentum equa- 
tion, however, we have neglected the collision termt£l in 
the limit t^n > (m e /mi)(i/ & i/cj)uii. We have further ne- 
glected the electron viscosity term in Eq.td). The re- 
maining equations are the energy equationsB, 
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\ n <^TT + ^( V II ' «*e||) = V,| • (x|V||T ) 



(9) 



•0.71V r (n e T e U||), (10) 



where Xi|| = 3.9nio7io/m;i'M and x e || = 3.2n e0 T e0 /r7i e ^ ei 
are the parallel ion and electron-thermal conductivities, 
respectively. In the above equations, the terms ~ u>~ e \ 
have neglecetd owing to the fact that ui < oJ C e,i- 

We consider a small electrostatic perturbation and lin- 
earize the above equations. The linearized equations are, 



dn el 
dt 



dug 
V||(n e0 M e ||i 



+ Vy(n 4O u i y 1 ) = 0, (11) 

n-elUejlo) = 0, (12) 



where equilibrium and perturbed quantities are defined 
by the subscripts and 1, respectively and u e || is the 
zeroth order electron-drift velocity with respect to the 



ions and dust grains. The electrons drift with respect to 
the ions and dust particles, so that Uio = Udo = 0. The 
ion motion along the magnetic field B is given by the 
linearized parallel momentum equation 



Out 



dt 



-V|| (enioVi 
- II Vtj Willi, 



riioTn + rinTio) 



(13) 



where tpi is the perturbed electrostatic potential. For 
lu < v e i, the parallel electron momentum equation be- 
comes 

= V|| (eneoPi - n el T e0 - n e0 T e i) ~ 0.71n e0 V\\T el 

- 0.51m e ^ el0 (n e0 M e ||i —Ucouqi) 

- 0.51m e n eO u e || O ^ e i, (14) 

where v e \ = I'eio (^ei /^eo — 3T e i/2T e o) and the equilib- 
rium drift velocity ti e M = —eEiiQ/m e v e i. The perturbed 
electron and ion temperature equations are given by 



3 dTn 2 



^n e0 ( + « e ||oV|| ) T e i + neoTeoVnUeH! = Xe||VjjT e l - 0. 71n e0 T e0 V|| (lt e || i - uqi) 



- 0.71u e || V|| (n e0 T el + n el T e0 ), 



(15) 



(16) 



Finally, we write the equations .for dust charge fluctua- 
tions as given by Jana et al.EZEa and the quasineutrality 
condition as 



d_ 

dt 

e(n e i - Tin) 



Vj Qdi 
ndoQdi 



n i0 



n e i 



0. 



(17) 
(18) 



where Qda = eZdo, Zdo is the equilibrium charge num- 
ber on the surface of the dust grains, Qdi = eZdi, Zdi 
is the charge fluctuation, n = e\I e \{T~^ ~ W ~ )/C with 
Wq = Tio — eipfo, (j)fo is the potential at the dust surface, 
|7 e | ~ Ii <~ ennyTta 2 c s is the equilibrium electron (or ion) 
current at the dust surface, and C ~ a being the grain 
capacitance. In writing Eq.(|l8|), we have neglected the 



dust density fluctuations as the ion transit tiaae (~ L/c s ) 
is much shorter than the dust transit timeE3 [~ L/cd*, 
Cd* ~ {Zj Q ndoTio/niomd) 1/2 ] along the magnetic field. 

It should, however, be noted that in writing the 
above equations, we have assumed a stationary back- 
ground equilibrium. On the other hand, a non-stationary 
background equilibrium may lead to stabilization of the 
ion-aco«stic instability in a low temperature collisional 
plasmaEl 

We now take perturbations of the form 



/(z,t)~/ ie -^- fc i z ) 



and write Eqs.(|l[|T|) as 



(19) 



3 Wjjo 
2 X 
3 



C||Mz||l 



fen Zt. 



e||l 



hum 



(u -W|| )' 

1,2 2 / 
ft || c s 

(w + i£||) 
.fen 



fl>i Ti 

<p + — + — 

T T 



1-71 +i--^-)T e -(p = i — {u e \\ x + u e \\ h e - iti||i), 

Xe 



+ iXi\\ T i = fc pi||l> 



(20) 
(21) 

(22) 

(23) 
(24) 
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— 2.21cj||q + ix e \\ I T e = fc||(1.71u e p — 0.71uj||i) + 0.71uj\\oh e 



(25) 



where r = T e0 /T i0 , h e = n e i/n e0 , hi = nn/n iQ , % = Tn/T i0 , T e = T el /T e0 , (p = e(pi/T e0 are the normalized variables 
with Xe — k 2 c 2 e /Q.51v eU u\\ = 0.96fc 2 c 2 /vu, uj\\ = fe||it e ||o, Xe\\ = 3.2/c 2 c 2 /V ei , and Xi\\ = 3.9fc 2 cf/^. We can now 
combine Eqs.(O) and (J18|) as 



l + i 



Zdon-do Iq 
n e0 (w + «??) 



+ 1- 



(26) 



where Jo = \Ie\/eZdo and Ad = (Teo/47rn e oe 2 ) 1 ^ 2 is the plasma Debye length. 

In the limit uj ~ U||q < Xe\\ i- e - 3.2(fc|| \ e ) 2 v e i/uj > 1, Eq.(p5|) can be re-written, using Eq.(pl|), as 



, (1.71^-^l| ) , w . 
-i — n e + zU.71 rii 

Xe\\ Xe\\ 



From Eqs.(|q) and (||), we get 



T, 



OJ 



Substituting Eqs.(||) and (||) into Eq.(|g), we obtain 

1.2 2 p 2 



Using Eqs.(^0|), (|2lJ), and fl27j), Eq.(||) can be written as 

' . w (1.71a; - oj[|o 
1 — « il.71 — 

Xe Xe\\ 



2„2 



kf,c 



n e -<p 



LU(UJ + i//||) 



-il.75— hi. 

Xe 



(27) 



(28) 



(29) 



(30) 



In the limit of ui > r), Eqs.(E6|), (129), and (BOh can be combined to obtain the following dispersion relation, 



w(w + ?;/i||) (w + iAn) + %n) 



l + i 



Zdon-do Iq 
n e0 lu 



»e0 



"(0 



2 , ,2 



n e0 w(w + «/i||) 



w n e0 



Xe 



il.75- — - ( l + i 

Xe n M 



For T, 



Z dan do Iq 
n e0 uj 

„2 n 2 



Xe\\ 

n l0 



(31) 



io ~ T e o and fcjjc 2 /w^ e i > 1, it is obvious that w > A|hXi|| (i- e - k^cj/uun < 1). We now write Eq.(|3l|) as 

e(/c|| , w) = and cj = u r + i"f with cj r > 7. Taking e = /t|| /cjf°' as an expansion parameter, the real part of u, to the 
first order, can be written as, 



0J r : 

where uj^ is given by 

The growth rate can now be written as 

7 - 



l.nkict 



2 2 ^11 n i0 W l|0 



uj1 0) n e0 Xe\\ 



(0 > - k«r 



n l0 5 T l0 
n-eo 3 T e o 



1/2 



(32) 



(33) 



2 9 w 2 

L.2 2 

n m «||C S 



- -7) 2- J 1 



1.71 



Wj|0 



Zdon d Q Xe\\ 
n i0 Xe 



(34) 
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FIG. 1. Normalized growth rate 7 of the ion-acoustic wave 
vs. normalized wave number k for different values of the pa- 
rameter S = n e o/nio- The curves labeled as (a), (&), and (c) 
are for 8=1 (no dust), 0.4, and 0.3, respectively. The other 



parameters are u = 6.0, I = 0, and A e 



0.01. 



FIG. 2. Behavior of 7 vs. k for (a) 7 = 0, (6) 0.01, (c) 0.05, 
and (d) 0.1. Dust density parameter 5 = 0.3 and A e ,i = 0.01. 



We note that, at the zeroth order, the phase veloc- 
ity of the ion-acoustic wave is modified with a factor of 

_i£ _|_ ^_l£ j xhc presence of dust (decreasing ra- 

tio of n e o/rijo) leads to an increase in the phase velocity 
of the ion-acoustic wave. 

In the expression for growth rate of the ion-acoustic 
wave Eq. (|34|) , the collisional damping is represented by 
the first two factors, the parallel ion-viscosity and parallel 
ion-thermal conductivity. The third term causes damp- 
ing due to charge fluctuation on the surface of the dust 
grains, which is proportional to the electron/ion current 
(Io) on the dust surface. We note that in absence of 
dust charge fluctuation (i.e. Iq — 0), increasing nega- 
tively charged dust density with respect to electron and 
ion density, results in an increase of the growth rate 



of the ion-acoustic instability. This can be seen from 
Fig.l, where we have plotted the normalized growth rate 
7 = jL/c g against normalized wave number k = k»L, as 
obtained from the solution of the full dispersion relation 
Eq. (|3l|) , at different negatively charged dust concentra- 
tion 5 = n e o/riio. Note that S = 0.3 means 70% of the 
electronic charges is now carried by the massive dust par- 
ticles and S = 1 means no negatively charged dust grains. 
The other parameters in Fig.l are u = it e ||o/c s = 6, 

A = Ai. e /L = 0.01, and / = I$L/c s = 0. As in ion- 
acoustic instability without dust particles, the positive 
growth rate, here also, appears only when the drifting 
electron velocity u e ||o? exceeds several times the phase 
velocity. Therefore, the necessary condition for required 
electron-drift relative to ions and negatively charged dust 
grains to excite the ion wave instability is 



u e \\o 



5 T i0 
3 T e o 



1/2 



(fc||Ae )2l^Zi£!^ 



1.8 



Zda n do \ 
nm J 

5 n i0 T e0 
3 



1 g7 / Z do n dQ \ n i0 /mo 5 T i0 

V n>i0 J n e0 V n e0 3 T e0 



-1 



■7Ta 2 \ e n l0 
Zdo 



1/2' 



(35) 



where we have used, only the zeroth order expression 
for Lu r . Note that the ion- viscosity and ion-thermal con- 
ductivity raise the critical electron current for the ion- 
acoustic wave to be unstable [the term in the square 
bracket in Eq.(|3^)], whereas the electron-thermal con- 
ductivity tends to increase the growth rate, as seen from 
Eq.©. 

The effect of increasing negatively charged dust den- 
sity on growth rate of the ion-acoustic instability may 
be compensated by the presence of charge fluctuation on 



dust surface, which has a damping effect, as shown in 
Fig. 2. In Fig. 2 we have taken the dust density parame- 
ter (5 = 0.3. The other parameters are same as in Fig.l. 
We show in Fig. 3, the behavior of 7 versus 8 for different 
values of dust current 1 = 0, 0.05, 0.1, and 1.0. Other 
parameters chosen are u = 6, A = 0.01, and k = 5.0. The 
overall behavior is consistent with Fig.l and Fig. 2. 
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k 

FIG. 3. Dependence of 7 on the dust density parameter S 
for different fluctuation level of the dust charge, (a) I — 
i.e. no dust charge fluctuation, (6) I = 0.05, (c) I = 0.1, and 
(d) I — 1.0. Note the change in behavior of 7 as the dust cur- 
rent increases. The normalized wave number k = 5.0, u = 6.0, 
and A e ,i = 0.01. 

However, it is important to note that, although a de- 
creasing 8 causes peaking up of the growth rate, at any 
given wave number, a rapid dust charge fluctuation rate 
(larger Jo) causes the wave to damp at higher dust con- 
centration, which is opposite in behavior to the case of 
absence of dust charge fluctuation (Fig.l and 2), as can 
be seen from all the curves in Fig. 3. 

For typical experimental parameters n e o/nio ~ 0.3, 
Zdon d o/n M ~ 0.7, T e0 - T i0 , and v ei /vu ~ (mj/m e ) 1/2 , 
the necessary threshold condition for instability reduces 
to 



u e||0 



> 2.24c< 



0.6 



1.7 

Tra 2 X e riio 



1/2 

— i (o.7/. t ;a- 



Z, 



d() 



(36) 



It is to be noted that the effects of charge fluctuation at 
the dust surface will tend to play an important role when 
the mean free paths are comparable to the plasma length 
along the magnetic field and the term 



1/2 



(37) 



Also notice that in high density and low temperature 
plasma limits, the effects ion-viscosity and thermal con- 
ductivity will be important, typically at wavelengths such 
that 



fciiA e ~ 1.5 



1/4 



(38) 



The condition ( p6|) can be realized in typical labora- 
tory situations^. For example, with an average size of 
dust grains to be few microns (a ~ 10 -4 cm) and a cor- 
responding Zdo ~ 10 4 , the condition (|36"|) yields 



*e||0 



> 4c s 



(39) 



for the excitation of the ion-acoustic instability, where 
we have taken the collisional parameters as \ e ,%/L ~ 0.01 
and k = k\\L ~ 2. The typical parallel electric field re- 
quired for the corresponding threshold electron current 
is - O.lV/m. 



IV. CONCLUSIONS 

We have studied the ion-acoustic instability in a col- 
lisional dusty plasma with fluid equations, where dust 
grains are treated as massive and negatively charged com- 
ponent in a multicomponent plasma. The regimes are 
clearly marked out where the theory is applicable, espe- 
cially in a relatively high density (10 10 cm -3 ) and low 
temperature plasma where the effect of collisions cannot 
be neglected. While treating the negatively charged dust 
particles, we take into account the effect of charge fluctu- 
ation on the dust surface. We have shown that in such a 
plasma, which are routinely produced in laboratory, there 
is a significant impact of electron-ion collisions, even at 
large mean free paths (i.e. A/ > L). We have derived the 
threshold electron-drift velocity required to drive the ion- 
acoustic instability. It is shown that the electron-thermal 
conductivity and the dust charge concentrations reduce 
the threshold value of electron current for driving the ion- 
acoustic mode. In particular, the ion-viscosity and ion- 
thermal conductivity raise the threshold current. And a 
similar effect due to dust charge fluctuations is also found 
when the mean free paths are of the order of the plasma 
length. 
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